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We use chiral perturbation theory for 5(7(2) to compute the leading corrections to the thermal 
mass of the pions and the pion decay constant in the presence of the magnetic field in a low- 
temperature expansion. The magnetic field gives rise to a splitting between M v o and M w ± as well 
as F n o and F 7T ± . We also calculate the free energy and the quark condensate to next-to- leading 
order. The results suggest that the critical temperature T c for the chiral transition is larger in the 
presence of a constant magnetic field, in agreement with most model calculations but in disagreement 
with recent lattice calculations. 



Introduction: Chiral perturbation theory (ChPT) pro- 
vides a systematic framework for calculating properties 
of QCD at low energies ChPT is not an expan- 

sion in powers of some small coupling constant, but it is 
a systematic expansion in powers of momenta p where 
a derivative counts as one power and the quark masses 
count as two powers. Chiral perturbation theory is a 
nonrenormalizable quantum field theory in the old sense 
of the word. This means that a calculation at a given 
order in momentum p, requires that one adds higher- 
order operators in order to cancel the divergences that 
arise in the calculations at that order. This implies that 
one needs more and more couplings and therefore more 
experiments to determine them. However, this poses 
no problem, as long as one is content with finite preci- 
sion. This is is the essence of effective field theory. The 
chiral Lagrangian that describes the (pseudo)Goldstone 
bosons are uniquely determined by the global symme- 
tries of QCD and the assumption of symmetry breaking. 
The Lagrangian C c s consists of a string of terms that in- 
volve an increasing number of derivatives or quark mass 
factors, each multiplied by a low-energy constant (LEC) 
li. However, QCD is a confining and strongly interacting 
theory at low energies. Thus the couplings U of the chi- 
ral Lagrangian cannot be calculated from QCD. Instead, 
the couplings are fixed by experiments 

The thermodynamics of a pion gas using ChPT was 
studied in detail in a series of papers 25 years ago 0- 
0]. The thermal pion mass and the thermal pion de- 
cay constant were evaluated at next-to-leading order, 
while the pressure and the temperature dependence of 
the quark condensate were calculated to next-to-next-to- 
leading order. In this letter, we present the calculation 
of these quantities in the presence of a constant mag- 
netic background. The details of the calculations can be 
found in Ref. Si. 



QCD in external magnetic fields has received a lot 
of attention in recent years due to its relevance in sev- 
eral physical situations. For example, large magnetic 
fields exist inside ordinary neutron stars as well as mag- 
netars [9|. In the latter case, the cores may be color 
superconducting and so it is important to study the ef- 
fects of external magnetic fields in this phase. Similarly, 
it has been suggested that strong magnetic fields are cre- 
ated in heavy-ion collisions at the Relativistic Heavy-Ion 
Collider (RHIC) and the Large Hadron Collider (LHC) 
and that these play an important role [l(|. In this case, 
the magnetic field strength has been estimated to be up 
to B ~ 10 19 Gauss, which corresponds to \qB\ ~ 6m^, 
where \q\ is the charge of the pion. This has spurred 
the interest in studying QCD in external fields. At zero 
baryon chemical potential this can be done from first 
principles using lattice simulations and some recent re- 
sult are found in [ljT - [l3j |. 

Chiral perturbation theory has been used to study the 
quark condensate in strong magnetic fields at zero tem- 
perature 14- 13], while the quark-hadron phase transi- 



tion was investigated in Ref. [18| . The effects of external 
magnetic fields on the chiral transition have been studied 
in detail using the NJL model 19|-|28|. the Polyakov- 
loop extended NJL model [1^, 3Q|, the quark-meson 
model 0, ED, IpK the (P)QM model [H 13, the 
linear sigma model [36[ and the MIT bag model [35j 

Chiral perturbation theory: As explained in the in- 
troduction, chiral perturbation theory is low-energy ef- 
fective field theory that can be used to systematically 
calculate physical quantities as a power series in mo- 
mentum. The effective Lagrangian is given by an infi- 
nite string of operators involving an increasing number 
of derivatives or quark masses. Schematically, we can 
write £ e ff = + + + ... where the superscript 
indicates the powers of momentum. 



By expanding the Lagrangian to fourth order in the pion fields 7Tj, we obtain 
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where we have defined the complex pion fields as tt± = ^(tti ± 1^2) and = BS^Xi- Similarly, expanding 
to second order in the pion fields yields 
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The Lagrangian is very complicated as it contains 
more than 50 terms for SU(2) [4\. However, only one 
term is relevant for the present problem (T5L nl\ , namely 
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We have used the parametrization U — e l7TiTi ^ F , where 
Tj are the Pauli matrices. This is different from the 
parametrization used in l]J 3] and so the expressions 
for C also differ. However, we get the same result for 
physical quantities as we should. 

The chiral Lagrangian comes with a number of un- 
determined parameters or low-energy constants. These 
parameters can be determined by experiments; however, 
loop corrections involve renormalization of them. The 
relation between the bare and renormalized parameters 
has been can be expressed as k = — [7 + 1 — h] , 

where 7, are coefficients and U are scale-independent pa- 
rameters 0. 

Pion mass and pion decay constant: The pion masses 
M^o and M T ± are defined by the position of the pole 
of the propagator. At leading order, the expressions for 
them are divergent and requires renormalization of the 
parameters l 3 , I5, and l&. The result is 
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where the pion mass M 2 in the vacuum is given by M 2 — 
MHl-j^wh), Ih=M 2 l og ML„ A f 2„ 2C (i,o )(0; i + 

x)\qB\, J 1= 4/3 2 f™ P 2 dp/^/p 2 + M 2 {ePVp 2 +M 2 - 1), 

Ji B =2EZ=ofo™ dp/V^TMl^V^l - 1), where 
M|> = M 2 + (2m + l)qB, and m denotes the Landau 
level. 



In order to calculate the pion decay constant, we 
need to evaluate the matrix elements (0|^l^|7r ) and 
(0|^4j l^^) ■ At zero magnetic field, these are identical, 
but there are two pion decay constants at finite tempera- 
ture; one for the time component and one for the spatial 
component of A u . The difference between them is an 



order-p effect 37]- The matrix elements are propor- 



tional to zp M and the prefactors are denoted by F^o and 
F„± , respectively. The expressions are divergent and 
requires renormalization of the Z4 and the renormalized 
result is 



F^o = F„ 

F„± = Fir 
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where the pion decay constant F^ in the vacuum is 
r T 4 ). Note that these are different in a 
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magnetic field since the contributions from the loops are 
different [8]. 

Free energy and quark condensate: We are inter- 
ested in the contributions to the free energy that are 
due to a nonzero magnetic field and finite tempera- 
ture. We therefore write the contribution to the free 
energy at the n'th loop order, T n , as a sum of three 
terms T n = J"™ c + + where F^ c is the free 
energy in the vacuum, i. e. B = T = 0, is the 
zero-temperature contribution due to a finite magnetic 
field, and J-^ is the finite-temperature contribution. The 
strategy is to isolate the term J-^ ac and subtract it from 
T n . This term contains ultraviolet divergences which 
are removed by renormalization of the low-energy con- 
stants of the chiral Lagrangian in the usual way and the 
renormalized .F™ represents vacuum energy of the the- 
ory. The sum + J 7 ^ generally contains ultraviolet 
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divergences as well and it is rendered finite by renormal- 
izing the Us. If we express the one-loop free energy in 



J 



terms of M„ instead of M, the dependence of I3 cancels 
in the expressions for Ti +2 - and J 7 ^^- This yields 
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where i T = M 2 /2|gB|, d(M 2 ) = 8(4:Tr) 4 d }i + ^(Iq — I5) log jjs, and A is the renormalization scale. We note that the 
NLO correction to the free energy in the chiral limit does not vanish in contrast to the case of zero magnetic field. 
At finite temperature, the quark condensate is (qq) = (0\qq\0)(l + -jh where the constant c is defined by 

c= -F 2 (0 1 qq | 0) ~ 1 7] . In the chiral limit , we have c = 1. In that case, the quark condensate reduces to 



(qq) = (0|gg|0Wl + 
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Results and discussion: We first notice that we in the 
limit B — > recover the temperature dependence of M^, 
Fn, J~, and (qq) as in p-lll- Similarly, we obtain the 
T = result for the free energy and the B dependence of 
the quark condensate as in [15l4l7j . We also notice that 
temperature dependence of the charged pion mass is the 
same as for vanishing magnetic field. The only difference 
is a temperature-independent constant proportional to 
(qB) 2 /F 2 . 

In the remainder we focus on the chiral limit. It is 
then easy to show that J\T 2 > jf\qB\ for all values 
of B and T. This implies that the pion decay con- 
stants F^o and F 7r ± are larger than F„. Moreover, for 
small values of | qB | , we can calculate the first corrections 
due to nonzero B as a power series in y/\qB\/T. One 
finds f> = F^l - T 2 /12F 2 + 5 y/\qB~\T / '48tt F 2 + ...) 
and F„o = F w ±(l - T 2 /12F 2 + 5^/\qB\T /96nF 2 + ...). 
Similarly, one can expand the quark condensate around 
\qB\ = and obtains (qq) = <0]gg|0)(l + ...T 2 /SF 2 + 
$yf\q~B\T/WiiF 2 + ...). 

In Fig. [TJ we show the quark condensate Eq. as 
a function of temperature for \qB\ = 900 (MeV) 2 (solid 



curve). For comparison, we also show the quark con- 
densate for \qB\ = (dashed curve). We are using the 
experimental value F n — 93 MeV. The quark conden- 
sate for vanishing B goes faster to zero than it does in 
the presence of a magnetic field. This indicates that the 
critical temperature T c for the chiral transition is higher 
in a nonzero magnetic field. Of course, this conclusion 
is cautious since the behavior of the quark condensate 
in the vicinity of T c is beyond the reach of the low- 
temperature expansion. The result is in line with some 
model calculations, both mean-field type 
and beyond [34| . This is in disagreement with the recent 
lattice simulations 13j where a significant decrease of 
the crossover temperature was observed. On the other 
hand, the lattice result in Ref. 11[ report a small in- 
crease of T c with B. Model calculations that seem to 
indicate a decrease can be found in Refs. 



18, 32, 3 



A possible explanation for the different behavior of T c 
as a function of B found in [32[ and [33| could be the 
different treatment of the fermionic vacuum fluctuations 
and related renormalization issues. 

In Ref. [l8| |. the authors use ChPT to investigate the 
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FIG. 1: Quark condensate in the chiral limit as a function 
of temperature for \qB\ = 900 (MeV) 2 (solid curve), and 
\qB\ = for comparison (dashed curve). 

quark-hadron phase transition as a function of the mag- 
netic field at the physical point. They compare the pres- 
sure in the hadronic phase with that of the quark-gluon 
plasma phase for an ideal gas of quarks and gluons, and 
subtracting the vacuum energy due to a nonzero gluon 
condensate {g 2 G flu G^ v ). For weak magnetic fields, the 
transition is first order. The line of first-order transi- 
tions ends at critical point. From this temperature on- 
wards, the transition is a crossover. Normally, however, 
the critical temperature is determined by the behavior 
of the quark condensate. At the physical point, the con- 
densate never vanishes and transition is a crossover. The 
crossover temperature is often defined by the inflection 
point of (qq) as a function of temperature. 

In conclusion, we have used chiral perturbation theory 
to calculate the pion masses, the decay constants, the 
free energy and the quark condensate at finite temper- 
ature in a magnetic background. Given the conflicting 
results for T c as a function of B of various model calcu- 
lations and lattice calculations, clearly more work needs 
to be done. 
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